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Abstract 

In this paper we study the cosmological evolution of the holographic dark energy in a cyclic uni- 
verse, generalizing the model of holographic dark energy proposed by Li. The holographic dark energy 
with c < 1 can realize a quintom behavior; namely, it evolves from a quintessence-like component 
to a phantom-like one. The holographic phantom energy density grows rapidly and dominates the 
late-time expanding phase, helping realize a cyclic universe scenario in which the high energy regime 
is modified by the effects of quantum gravity, causing a turnaround (and a bounce) of the universe. 
The dynamical evolution of holographic dark energy in the regimes of low energy and high energy is 
governed by two differential equations, respectively. It is of importance to link the two regimes for 
this scenario. We propose a link condition giving rise to a complete picture of holographic evolution 
of a cyclic universe. 
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The astronomical observations over the past decade imply that our universe is currently dominated 
by dark energy, which leads to an accelerated expansion of the universe (see e.g. Refs. [U [21 13]). The 
combined analysis of cosmological observations suggests that the universe is spatially flat and consists 
of about 70% dark energy, 30% dust matter (cold dark matter plus baryons), and negligible radiation. 
The basic characteristic of dark energy is that its equation of state parameter w (the definition of w 
is w = p/p, where p is the energy density and p is the pressure) has a negative value [w < —1/3). 
The most obvious candidate for dark energy is the cosmological constant A [4] for which w = — 1 (for 
reviews see e.g. Refs. [S]). However, the cosmological constant always suffers from the "fine-tuning" 
and "cosmic coincidence" problems. Another candidate for dark energy is the energy density associated 
with a dynamical scalar field; a slowly varying, spatially homogeneous component. Typical examples of 
such a type of dark energy are the so-called "quintessence" [6] and "phantom" [7] models. Quintessence 
dark energy provided by a canonical scalar field has an equation of state w > — 1, while the phantom 
energy associated with a scalar field with a negative-kinetic energy has an equation of state w < — 1. It 
is remarkable that for phantom dark energy in this scenario all the energy conditions in general relativity 
(including the weak energy condition) are violated. Due to such a supernegative equation of state, the 
phantom component leads to a "big rip" singularity at a finite future time, at which all bound objects 
will be torn apart. 

In the phantom scenario, generically, there exist two space-time singularities in the universe; one is the 
initial "big bang" singularity, the other is the future "big rip" singularity. The space-time singularities are 
disgusting for the majority of theorists; thus a mechanism for avoiding the initial and future singularities 
are attractive for physicists and cosmologists. An effective way for eliminating the singularities is to 
introduce a p 2 term with a negative sign to the Friedmann equation if the energy is very high. Such 
a modified Friedmann equation with a phantom energy component leads to a cyclic universe scenario 
in which the universe oscillates through a series of expansions and contractions. Phantom energy can 
dominate the universe today and drive the current cosmic acceleration. Then, as the universe expands, it 
becomes more and more dominant and its energy density becomes very high. When the phantom energy 
density reaches a critical value, a very high energy density, the universe reaches a state of maximum 
expansion, which we call the "turnaround point" , and then it begins to recollapse, according to the 
modified Friedmann equation. The contraction of the universe makes the phantom energy density dilute 
away and the matter density dominate. Once the universe reaches its smallest extent, the matter density 
hits the value of the critical density, the modified Friedmann equation leads to a "bounce" , making the 
universe once again begin to expand. 

The idea of an oscillating universe was first proposed by Tolman in the 1930's [8]. In recent years, 
Stcinhardt, Turok and collaborators [3] proposed a cyclic model of the universe as an alternative to the 
inflation scenario, in which the cyclicity of the universe is realized in the light of two separated branes. 
The cyclic scenario discussed in this paper is distinguished from the Steinhardt- Turok cyclic scenarios in 
that the phantom energy plays a crucial role. In the oscillating or cyclic models, the principal obstacles 
against success come from the problems of black holes and entropy. For the discussions of the problems 
of black holes and entropy in the cyclic scenario, see e.g. Refs. [TUlfTT] . 

Usually, the phantom energy density becomes infinite in a finite time, leading to the big rip singularity. 
However, we expect that an epoch of quantum gravity sets in before the energy density reaches infinity. 
Therefore, we arrive at the notion that quantum gravity governs the behavior of the universe both at 
the beginning and at the end of the expanding universe, where the energy density is enormously high. 
This high energy density physics may lead to modifications to the Friedmann equation, such as in loop 
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quantum cosmology |12j and braneworld scenarios |13[ 114] , which causes the universe to bounce when it is 
small, and to turn around when it is large. At high energy densities, we employ the modified Friedmann 
equation 



where H = a/a is the Hubble parameter, G is the Newton gravity constant, and p c is the critical 
energy density set by quantum gravity, distinguished from the usual critical density 3M^H 2 (where 



the effective theory of loop quantum cosmology [12] , and also from the braneworld scenario [2] . In loop 
quantum cosmology, the critical energy density can be evaluated as p c w 0.82p p i, where p p \ = G~ 2 = 
2.22 x 10 76 GeV is the Planck density. In the braneworld scenario, p c = 2a, where a is the brane tension, 
and a negative sign in Eq. (JT|) can arise from a second time-like dimension; but that gives difficulties 
with closed time-like paths. In models motivated by the Randall-Sundrum scenario [13] . the most natural 
energy scale of the brane tension is of the order of the Planck mass, but the problem can generally be 
treated for any value of a > TeV 4 . Once the energy density of the universe reaches the critical density 
p c , the universe changes its evolution direction. At that energy scale, if it has been expanding, it turns 
around and begins to contract; if it has been contracting, it bounces and begins to expand. Modifications 
to the Friedmann equation thus motivate the bounce and the turnaround, both of which are nonsingular. 

Recently, considerable interest has been stimulated in explaining the observed dark energy by the 
holographic dark energy model. The holographic dark energy model is an attempt to probe the nature 
of dark energy within the framework of a fundamental theory originating from some considerations of 
the features of quantum gravity theory. Concretely speaking, this model is constructed in the light of 
the holographic principle [15] of quantum gravity. In the holographic scenario, the dark energy is a 
dynamically evolving vacuum energy density that can realize the phantom behavior. If the holographic 
dark energy becomes a phantom, this scenario will involve a big rip singularity in the far future unless 
the Friedmann equation gets a quantum gravity correction in the high energy regime as shown in Eq. (fT]). 
Thus, a phantom-like holographic dark energy can play a crucial role in realizing a cyclic universe scenario, 
and inversely, a cyclic universe can endow the holographic dark energy with peculiar features. In this 
paper, we shall study the cyclic universe with a holographic phantom and investigate the characteristic 
of the holographic dark energy in such a cyclic universe. 

According to the holographic principle, the number of degrees of freedom for a system within a finite 
region should be finite and should be bounded roughly by the area of its boundary. In the cosmological 
context, the holographic principle will set an upper bound on the entropy of the universe. Motivated by 
the Bekenstein entropy bound, it seems plausible that one may require that for an effective quantum field 
theory in a box of size L with UV cutoff A, the total entropy should satisfy S = L 3 A 3 < Sbh = ^M^L 2 , 
where Sbh is the entropy of a black hole with the same size L. However, Cohen et al. [16] pointed out 
that to saturate this inequality some states with Schwartzschild radius much larger than the box size 
have to be counted in. As a result, a more restrictive bound, the energy bound, has been proposed to 
constrain the degrees of freedom of the system, requiring the total energy of a system with size L not to 
exceed the mass of a black hole with the same size, namely, L 3 A 4 = L 3 pa < LM 2 V This means that the 
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maximum entropy is in the order of Sg H . When we take the whole universe into account, the vacuum 
energy related to this holographic principle is viewed as dark energy, usually dubbed holographic dark 
energy. The largest IR cut-off L is chosen by saturating the inequality, so that we get the holographic 




(1) 



M p i = l/V 87rG is the reduced Planck mass). This modified Friedmann equation can be derived from 
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dark energy density 

p A = 3c 2 M^L- 2 , (2) 

where c is a numerical constant (note that c > is assumed), and as usual M p \ is the reduced Planck 
mass. Hereafter, we will use the unit M v \ = 1 for convenience. It has been conjectured by Li [T7] that 
the IR cutoff L should be given by the future event horizon of the universe, 

. , 7 dt' 7 da' 

RM = a J W)= a J ( ) 

t a 

Such a holographic dark energy looks reasonable, since it may provide simultaneously natural solutions 
to both dark energy problems, as demonstrated in Ref. |17j . The holographic dark energy model has 
been tested and constrained by various astronomical observations [18, 19J. For other extensive studies on 
the holographic dark energy, see e.g. Refs. [20| [21] . 

The holographic dark energy scenario reveals the dynamical nature of the vacuum energy. When 
taking the holographic principle into account, the vacuum energy density will evolve dynamically. The 
dimensionless parameter c plays a crucial role in the holographic evolution of the universe. As has been 
pointed out in Refs. [TH], the value of c determines the destiny of the holographic universe. When c > 1, 
the equation of state of dark energy will evolve in the region of —1 < w < —1/3. In particular, if 
c = 1 is chosen, the behavior of the holographic dark energy will be more and more like a cosmological 
constant with the expansion of the universe, such that ultimately the universe will enter the de Sitter 
phase in the far future. When c < 1, the holographic dark energy will exhibit a quintom-likc evolution 
behavior (for "quintom" dark energy, see, e.g., Refs. [22] and references therein), i.e., the holographic 
evolution will make the equation of state cross w = — 1 (from w > — 1 evolves to w < —1). That is 
to say, c < 1 makes the holographic dark energy today behave as a phantom energy that leads to a 
cosmic doomsday ("big rip") in the future. Nevertheless, as discussed above, at high energy densities the 
Friedmann equation may be modified to Eq. ([T]) due to some possible quantum gravity effects, which 
can successfully eliminate the big rip singularity. It is remarkable that the analyses of the observational 
data imply that the value of c in the model of holographic dark energy is very likely less than 1 [TB], i.e., 
the holographic dark energy is very possibly behaving as a phantom energy presently. Intriguingly, then, 
considering the modified Friedmann equation |T]) in the high energy regime, the holographic dark energy 
(with c < 1) along with some dust-like matter components can realize a cyclic universe scenario in which 
the cosmological evolution is nonsing ularQ 

First, consider the low energy regime of the universe, p <C p c - At this regime, for the universe we 
have the usual Friedmann-Robertson- Walker (FRW) case, 3H 2 = p. Consider a universe filled with a 
matter component p m (including both baryon matter and cold dark matter) and a holographic dark 
energy component pa; then the Friedmann equation reads 

3H 2 = Pnl + p A . (4) 



*In general, cyclic universe models confront two severe problems making the infinite cyclicity impossible. First, the 
black holes in the universe, which cannot disappear due to the Hawking area theorems, grow ever larger during subsequent 
cycles, and they eventually will occupy the entire horizon volume during the contracting phase, so that the calculations 
break down. The second problem is that the entropy of the universe increases from cycle to cycle due to the second law 
of thermodynamics, so that extrapolation into the past will lead back to an initial singularity. In this paper, we do not 
consider these difficulties. 
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Defining the fractional densities £l A = p\/3H 2 and O m = p m /3H 2 = Q^HqH 2 a 3 , where a is the scale 
factor of the universe and ao = 1 has been set, the Friedmann equation (j4|) can also be rewritten as 

^ = n° m a- 3 + n A ^. (5) 

Combining the definition of the holographic dark energy @ and the definition of the future event horizon 
we derive 

d In a c 

(6) 



Ha' HaJVl^ 
We notice that the Friedmann equation ([5|) implies 

— = Jail - fil) 1 -=, (7) 

Ha V 1 'HoJS^ 

Substituting into I©, one obtains the following equation: 



f°° e x '/ 2 JT^dx' = ce x ' 2 (8) 

where x = In a. Then taking the derivative with respect to x in both sides of the above relation, we easily 
get the dynamics satisfied by the dark energy, namely the differential equation for the fractional density 
of the dark energy, 

n' A = n A (i-n A ) (i + -^nf) , (9) 



c 

where the prime denotes the derivative with respect to x — In a. This equation describes the behavior of 
the holographic dark energy in the low energy regime completely, and it can be solved exactly |17j : 

In fl A - ln(l - JtT a ) + ln(l + y^I) - -t-^-z He + 1JtT a ) = In a + y Q , (10) 

z + c 2 — c 4 — c z 

where yo is an integration constant that can be determined by setting today as an initial condition, 



y = ln(l - n° m ) - ^- ln(l - JT^) + ^ + V^l) - Hp + 2y/T^). (n) 

From the energy conservation equation of the dark energy, the equation of state of the dark energy can 
be given [T7] : 

1 din OA 1 . 2 rrr . 

w = -l- -—£± = - -1 + -Vol . (12) 
3 a ma 3 c 

Note that the formula pa = [fl A /(l ~ fl A )]p® n a~ 3 and the differential equation of il A , see ((9]), are used in 
the second equality sign. 

As time passes by, the dark energy gradually dominates the evolution of the universe, Q A increases to 
1, and the most important term on the left-hand side of (flTJ|) is the second term; thus for large a, we get 

I 2 + c , s 8 2 + c 2 + c 

VUi = l-2-(2 + c)^M)e^-V-. (13) 
Since the universe is dominated by the dark energy for large a, we have 

Thus, using Eq. (|13|) in the above relation, we derive 

2(2 + c) 8 2 + c n 2(l-c) 

p A = 2-T^(2 + c)"^^)e— "Vm^ - ^ • (15) 
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It can be clearly seen from the above expression that the value of c plays a significant role in determining 
the hnal evolution of the dark energy. When c = 1, the holographic dark energy will become a cosmological 
constant that is related to p® n through the above relation without a. The choice of c > 1 makes the density 
of the dark energy continuously decrease, just like the case of quintessence dark energy. When c < 1, the 
density of dark energy ceaselessly increases with the expansion of the universe. The phantom behavior 
(c < 1) results in the density of holographic dark energy becoming enormously high when a is very large. 
As the universe goes into the high energy regime, quantum gravity begins to operate, giving rise to the 
modified Friedmann equation, JTJ). However, it is very hard to justify the borderline between the usual 
Friedmann equation and the modified Friedmann equation. Namely, one cannot accurately say when the 
usual Friedmann equation is replaced by the modified one. Hence, we have to set a criterion by hand. One 
can assume that quantum gravity begins to operate when pa = VPc, where, say, r] ~ 0(1CP 3 ) — 0(1O -2 ). 
We thus derive, at that moment, the scale factor of the universe, 



r m 



772^=^(2 + c) =(<=- 2 >e 



2(1- = ) 

(16) 



We assume that when a > a v the evolution of the universe is governed by Eq. (T]). Here, any other 
forms of energy have already been decayed away, and the dark energy is thus the unique component of 
the universe; so we have 

3if 2 =p A (l-^). (17) 

From this relation, we have 

n A = i + — — . (is) 

Pc - PA 

This indicates that fl\ > 1 when a > a v , even though the space of the universe is assumed to be flat. 
In concrete terms, when p\ <C p c , we have J7a - * 1 + ; when p\ — > p Cl we have Q\ — > 00. And, as a 
contrast, see Eq. (|13[) . for large a but a < a v , we have J7a — > 1~ . Now the question naturally arises of 
how to realize the tt\ = 1 crossing. This question is actually equivalent to the one of asking when the 
usual Friedmann equation should be replaced by the modified one. The transition of the two phases is 
ambiguous so that we have to set a connection when p\ = r]p c . Hence, the initial stage for Eq. (117|) is 
from Q\ = 1 + e, where e is a small positive number. 
The modified Friedmann equation can be rewritten as 

= * = ^1 (19) 

where the dimensionless parameter h is positive for an expanding universe and negative for a contracting 
universe. Here, since an expanding universe is considered for illustration, we take the positive value to 
h. Combining the definitions of holographic dark energy and future event horizon, namely Eqs. (JSJ) and 
©, yields 

' ,hl "' ' (20) 



Following Eq. (fT§|) . we have 
Substituting (21]) into (20|) yields 



ha' hay/Q\ 
1 v^A 



ha a^/flA — 1 



(21) 



dx' , " A = CVIIX , (22) 
e-VOv - 1 eVO\ - 1 
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where x = In a. Then, taking the derivative with respect to x in both sides of this equation, one obtains 
a differential equation for the fractional density of the holographic dark energy, 

n' A = 2fi A (Q A - 1) Q ^ - ij , (23) 

where the prime denotes the derivative with respect to In a. This differential equation governs the 
holographic evolution of the universe for the high energy regime. Note that here c < 1, and CI a > 1; 
hence f2 A is always positive, namely the fractional density of dark energy increases in time, the correct 
behavior as we expect. This equation can be solved exactly, and the solution is 

]ny/Shi+ C - ln(V^A-l)- 57^ r HV^A + 1) ~ ~, « ln(V»I - c) = In a + y v , (24) 

2(1 — c) 2(1 + c) 1 — c z 

where y v is an integration constant that can be determined by an appropriate initial condition. Now let us 
deduce the equation of state for the holographic dark energy in the high energy regime (a > a v ). Following 
the energy conservation equation of dark energy, we have w = —1 — (l/3)(dhip\/dhia). Writing 

Pa = 3flAh 2 p c = — ^- — -p c , (25) 
' 'a 

one can easily obtain 

w = ~(l+*y/Slf\. (26) 

Interestingly, this relation is the same as in the low energy regime; see Eq. (fT"2"|) . Note that here CI a is 
governed by Eq. (l2"3"j) . 

It has been pointed out that the borderline between the usual Friedmann equation and the modified 
one is rather ambiguous. One has to designate a contrived link condition; for example, one can assume 
that the transition has occurred when p\ = r/p c , where r\ ~ O(10~ 3 ) — 0(1O -2 ), say. Therefore, the initial 
stage of the high energy regime is contrived to be from a = a v , where ft\ = 1 + e with e = 77/(1 — r/). At 
the moment of a = a n and Qa — > 1 + , the most important term on the left-hand side of (|24| is the second 
term; we thus can determine the integration constant: 

y « = W^cj ln W^) - 2(TT^) ln2 - ih^ ln(1 - c) - lna "' (27) 

where a v is given by Eq. (fT6|) . It is remarkable that for a large 17a, when the universe approaches the 
turnaround point, the left-hand side of (|2"4"|) goes to zero. This gives the scale factor corresponding to the 
turnaround, 

a max = e-"', (28) 

where y n is given by Eq. (|27l) . Namely, the maximum scale factor (at the turnaround) of the universe 
in the cyclic universe scenario with holographic dark energy is totally determined by the constant y n . It 
should be noted that the link condition pa = rjp c is rather important in this scenario, even though it is 
somewhat artificial. The reason for contriving such a link condition is that the borderline between the 
classical and quantum gravities is not so clear. Next, let us give several numerical examples. We take 
c = 0.8, f^ n = 0.27, and h — 0.72 (here h is the dimensionless Hubble parameter of today), which gives 
rise to yo = —1-54. The critical density of the universe depends on the theory we use; if we use an effective 
theory of loop quantum cosmology, we have p c « 0.82p p i = 1.82 x 10 76 GeV ; if we use a braneworld 
scenario, we can treat the value of p c from TeV 4 to 10 76 GeV . Here we take the loop quantum cosmology 
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for illustration. The present density of dust matter is = 1.13 x 10~ 47 GcV 4 . Then, we can determine 
the values of and a max if the value of r\ is given. We only show two examples, r\ = 10 -3 and r; = 1CP 2 . 
The choice of rj = 10~ 3 gives a n = 2.86 x 10 240 and a max = 1.53 x 10 245 ; the choice of rj = 10~ 2 results 
in a„ = 2.86 x 10 242 and a max = 1.50 x 10 245 . 

In fact, this scenario has a fatal flaw, namely, a cyclic universe has no a future event horizon in 
principle, since an observer can eventually see the whole universe due to the cyclicity of the universe, if 
he/she waits a sufficiently long time. Therefore, the calculations in this paper break down in this regard. 
However, we can rescue the model by reconsidering the IR cutoff of the universe. In the original work of 
the holographic dark energy [T7], the choice of the future event horizon as an IR cutoff is only a conjecture 
for ensuring the acceleration of the universe. Now that the future event horizon cannot exist in a cyclic 
universe, we might as well make a modification to the future event horizon. We can define a "finite-future 
event horizon" , by replacing the infinity with a time T in the upper limit of integration in Eq. ([3]) , where 
T denotes the time of the turnaround. Choosing the finite-future event horizon as the IR cutoff of the 
universe undoubtedly makes the holographic dark energy meaningful in the cyclic universe, at least in 
the expanding stage. However, for the contracting stage of the cyclic universe, the IR cutoff is ambiguous 
for us. We can steer clear of this difficulty by considering the contracting stage as a time-reversal course 
of the expanding stage. It should be admitted that this assumption is rather strong. It should also be 
noted that the emphasis of this paper is placed on the holographic evolution in the high energy regime 
in the expanding branch and on how to link the low and high energy regimes. 

When the holographic phantom density reaches the critical density p c , the universe starts to turn 
around and contract. In the contraction phase, the physical rules are assumed to be totally the same as 
in the expansion phase, i.e., the high energy regime is governed by Eq. (|17[) and the low energy regime is 
dictated by Eq. (jl]). The only difference is that the Hubble parameter has a negative value, but this does 
not affect the evolutionary rules of the holographic dark energy; i.e., the dynamical evolution of the holo- 
graphic dark energy is still controlled by the differential equations © and (j2"3"|) . As the universe contracts, 
at first the energy density of the universe decreases, because the holographic phantom density decreases 
in importance (then, the quintom behavior makes the holographic dark energy become a quintessence-like 
component whose density increases as the universe contracts, although not as strongly as matter and ra- 
diation do), but then it again increases as matter and radiation become dominant. Eventually, the energy 
density reaches the high values at which the modifications to the Friedmann equation become important. 
Once the energy density again hits the same critical density p Cl the universe stops contracting, bounces, 
and once again expands. The bounce looks like a "big bang" for us, and at this point the universe has 
its smallest extent (smallest scale factor a m i n ) and largest energy density (the critical density p c ). An 
inflationary period may occur if the inflaton field can be excited by some mechanism that can solve the 
flatness and horizon problems, etc., and this can also generate the scale-invariant primordial perturba- 
tions seeding the structure formation. For the inflationary universe in a loop quantum cosmology, see, 
e.g., Refs. 23J. As the universe expands, its density deceases, and it goes through the radiation domi- 
nated and matter dominated periods, with the usual primordial nucleosynthesis, microwave background, 
and large structure formation. Around a redshift z ~ 0(1), the universe begins to accelerate due to 
the existence of dark energy (the holographic dark energy in this paper). The holographic dark energy 
with c < 1 can help realize the turnaround discussed above. It is noteworthy that the cyclic universe 
discussed in this paper is an ideal case, and there are still several severe obstacles existing in the cyclic 
cosmology, such as the density fluctuation growth in the contraction phase, black hole formation, and 
entropy increase, which can obstruct the realization of a truly cyclic cosmology. These problems are not 
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addressed in this paper. 

To summarize, in this paper we investigated the holographic dark energy in a cyclic universe. We 
generalized the model of holographic dark energy proposed in Ref. 17] to the case of a cyclic universe, 
and we studied the cosmological evolution of the holographic dark energy in such a universe in detail. The 
holographic dark energy with c < 1 can realize a quintom behavior; namely, it evolves from a quintessence- 
like component to a phantom- like component. The phantom energy density will become very large in 
the far future, which leads to a "big rip" singularity at a hnitc time in which all bounded objects are 
finally disrupted. However, when the energy scale becomes enormously large, quantum gravity effects 
may bring significant modifications to the Friedmann equation, leading to the possibility of the avoidance 
of singularities. A modified Friedmann equation, ([T]), along with a phantom component and other matter 
components can realize a cyclic universe scenario in which the cosmic evolution is nonsingular. In such 
a cyclic scenario, the large densities cause the universe to bounce when it is small, and to turn around 
when it is large. We investigated the holographic phantom cosmological evolution (c < 1) in such a cyclic 
universe in detail. The dynamical evolution of holographic dark energy in the low energy regime and in 
the high energy regime are rather different; they are governed by two differential equations, respectively. 
Linking the two regimes together is a very important mission for this scenario. We proposed a link 
condition connecting the regimes of low energy and high energy, which gives rise to a complete picture 
of the holographic evolution of the cyclic universe. 
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